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Quantum walks exhibit many unique characteristics compared to classical random walks. In the classical 
setting, self-avoiding random walks have been studied as a variation on the usual classical random walk. 
Here the walker has memory of its previous locations and preferentially avoids stepping back to locations 
where it has previously resided. Classical self- avoiding random walks have found numerous algorithmic 
applications, most notably in the modelling of protein folding. We consider the analogous problem in the 
quantum setting - a quantum walk in one dimension with tunable levels of self- avoidance. We complement 
a quantum walk with a memory register that records where the walker has previously resided. The walker is 
then able to avoid returning back to previously visited sites or apply more general memory conditioned 
operations to control the walk. We characterise this walk by examining the variance of the walker's 
distribution against time, the standard metric for quantifying how quantum or classical a walk is. We 
parameterise the strength of the memory recording and the strength of the memory back-action on the 
walker, and investigate their effect on the dynamics of the walk. We find that by manipulating these 
parameters, which dictate the degree of self- avoidance, the walk can be made to reproduce ideal quantum or 
classical random walk statistics, or a plethora of more elaborate diffusive phenomena. In some parameter 
regimes we observe a close correspondence between classical self- avoiding random walks and the quantum 
self-avoiding walk. 
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Quantum walks 1 " 4 , the quantum equivalent of classical random walks, have been studied extensively for 
their applications in quantum information processing 5 . Here a walker (a particle such as a photon) 
resides at a vertex (e.g. an optical mode) in a graph and is allowed to 'hop' along the edges in the graph to 
reach other vertices. In the classical random walk this process takes place randomly, whereas in the quantum case 
the walker coherently enters a superposition across different vertices. Numerous elementary optical demonstra- 
tions of quantum walks have been performed 6 " 14 , experimentally confirming the unique behaviour of quantum 
walks compared to classical random walks, and providing an alternate route towards optical quantum informa- 
tion processing 1516 . 

In this paper we consider one-dimensional quantum walks, where the walker has memory of its previous 
location history and the coin operator is a function of the memory. This can be used, for example, to preferentially 
avoid previously visited vertices or implement more elaborate conditional coins, allowing for tunable degrees of 
self- avoidance in the evolution of the walk. Self- avoiding walks have been studied extensively in the classical 
case 17 " 22 , having been applied to applications such as protein folding and percolation theory. 

Previous authors have considered quantum walks complemented by memory, where the memory is of previous 
coin values 23 " 26 . Here, however, we consider the case where the memory is of previous positions rather than of 
previous coins. 

Self-avoiding quantum walks may find applications in situations where ordinary quantum walks have been 
applied, but where we explicitly wish to avoid re-exploring previously visited sites. While classical self-avoiding 
walks have been applied to modelling protein folding 27 , long-chain polymers 28 and molecular conformation, 
perhaps quantum self- avoiding walks may be applied to modelling the analogous problems in the quantum 
context. 

In our model each vertex in the graph is complimented by a qubit, which marks whether the site has previously 
been visited. As the walker walks it leaves a record on the site at which it is presently located. The coin, which 
determines the subsequent step direction, is biased to avoid walking onto previously marked sites. We will see that 
due to the unitarity of our model, we cannot maintain a permanent record of where the walker has been and the 
walker has some chance of revisiting previously visited sites. 
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We find that self- avoiding quantum walks exhibit diverse diffusive 
characteristics, ranging from ideal classical to ideal quantum walk 
dynamics, and a plethora of more elaborate behaviours. In particular, 
we find that there is good agreement between the partially avoiding 
quantum walk and a corresponding partially avoiding classical walk 
in one-dimension. In the limit of maximum self- avoidance, the walk 
can exhibit significantly faster rate of spread than an ordinary 
quantum walk. 

Classical self-avoiding random walks. Self-avoiding classical ran- 
dom walks 17 " 20 are a variation on the standard classical walk, where 
the walker is unable to revisit any sites it has previously visited. A one 
dimensional, classical, completely self- avoiding walk is trivial; once 
the first step has been taken by the walker, every step thereafter has 
only one possibility. The variance of the walker will always be zero, as 
after the initial step there is only one possible location for the walker 
to be at any given step. Self-avoiding random walks become less 
trivial and more useful in two and three dimensions, but they also 
become exceedingly difficult to analyse in higher dimensions. 

Self-avoiding walks have the property that they explore a graph 
more rapidly, since evolution time is not wasted re-exploring already 
visited sites. Self- avoiding classical random walks are particularly 
useful in two and three dimensions as they find applications in mod- 
elling solvents and polymers 28 , protein folding 27 , molecular con- 
formation, percolation lattices, and drug discovery. When using 
classical random walks to model the proliferation of forest fires or 
the propagation of liquids through a porous material - two arche- 
typal applications for classical random walks - we wish to determine 
whether a route through a graph exists. Thus, there is no benefit in 
exploring vertices that have previously been explored and a self- 
avoiding walk will be far more efficient. 

Partially self-avoiding classical walks in one dimension have been 
studied previously and we review a standard model for one of these in 
Section below and find good comparison with the partial avoiding 
quantum walk. 

Quantum walks. A (discrete-time) quantum walk is a bipartite sys- 
tem comprising position and coin degrees of freedom, with Hilbert 
space H = / Hp®TLc- The position parameter specifies the location 
(vertex) of the walker in the graph, whilst the coin parameter specifies 
the direction the walker is heading. On a one- dimensional graph the 
state of the quantum walker at time t is of the form \\j/(t)) = Z Xj c ot Xy c 
(t) \x, c), where xeZ is the discrete position of the walker, and c = ± 1 
is the direction of the walker (left or right respectively). 

The evolution of the walk is specified by unitary coin (C) and 
step (s) operators, which have the actions U C j\xj) 

and \x,c)\-^>\x + c 9 c) respectively, where U is the unitary coin matrix. 
Thus the coin operator coherently manipulates the direction of the 
walker, leaving the position unchanged, whilst the step operator uses 
the coin parameter to update the location of the walker, leaving the 
coin unchanged. The evolution of the walk over t time steps proceeds 

as|^(0) = (S-C)V(0)). 

One of the interesting features of the quantum walk is that it 
exhibits a quadratically faster rate of spread (measured using vari- 
ance) compared to the classical random walk. The variance is defined 
as, o 2 (t) = H x p x (t) ' (x — fi(t)) 2 , fi(t) = H x p x (f) ' x, where p x (t) is the 
total probability of finding the walker at location x at time U given by 
summing over the coin degree of freedom, p x (t) = 2 c |a X; c {t)\ 2 , 
effectively tracing out the coin. 

Quantum self-avoiding walks. We fashion our self-avoiding quan- 
tum walk with three components: (1) we record the current position 
of the walker into the local quantum memory. This is performed 
using a 'memory update operator'; (2) we perform the usual coin 
operator, where the bias of the coin is determined by the state of the 
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Figure 1 | Illustration of the different subsystems in the self-avoiding 
quantum walk. The walker resides on a linear graph structure with a two- 
level memory qubit associated with each vertex, indicating whether the site 
has been previously visited. 

memory registers at the current location of the walker; (3) the usual 
step operator is applied, which results in the walker's movement 
being a function of whether the site has previously been visited. 
This allows, for example, the walker to reverse direction upon reach- 
ing a previously visited site. An illustration of the components in the 
self-avoiding quantum walk are shown in Fig. 1. 

Recording the current position — To model the self-avoiding 
quantum walk whereby the walker has memory of its previous loca- 
tion history and avoids positions where it has previously resided, we 
complement the walker by a sequence of qubits - one for each site in 
the lattice - where the qubits indicate whether the walker has prev- 
iously visited the respective site. In this case the state of the system at 
time t is of the form, 

l<K f )>= ^,c,q w --,to(t)\x,c,qi,...,q N ), (1) 

x,c,qi,-,q N 

where q { = {0, 1} is the qubit associated with site i, and there are N 
sites. Initially q l ■ = 0 |i i, indicating that none of the sites have been 
visited. Evidently, the size of the Hilbert space grows exponentially 
with the size of the lattice (\H\ =N-2 N+1 ), in contrast to linear 
growth for an ordinary quantum walk (\H\ =2N). This limits our 
numerical investigations to relatively small evolution times. 

We introduce a memory update operator, which updates the 
memory qubits as a function of the walker's position. We have cho- 
sen to implement the memory update operator unitarily, which 
implies that the memory is reversible and temporary, and that the 
record will evaporate after a second visit to the respective site. A 
permanent memory recording would require a Lindblad model, 
which would introduce decoherence into the system. Taking our 
memory recording to be unitary we fix it to be of the form, 

M(0 M )\x,c,qi, . . . ,q N )\-^R {x) (6 M )\x,c,q u . . . ,q N ), (2) 

a Pauli X rotation applied to the xth qubit. Now 9 M determines the 
strength of the action of the position on the memory. In the extreme 
case where 0 M = n/2, if the walker has not previously visited site x, 
the memory update operator flips the respective qubit to mark that 
this position has now been visited. Of course, a second visit to the site 
will flip the qubit again, marking that the site has not been visited. 
Note that this latter property is unavoidable if the evolution is to 
remain unitary. If 6m = 0 the position has no effect on the memory 
and the memory only couples with position via the coin operator. 

Biasing the coin according to neighbouring memory registers — The 
coin operator is defined analogously to the usual quantum walk coin 
operator, 

C(6c,0B)\x,c,qi, . . . ,<2at)i-> 

J2 U^\xlq u ...,q N ), (3) 
;=±i 

where is a 2 X 2 coin matrix that depends on the state of the 
qubit at location x, q x . This operator coherently manipulates the coin 
degree of freedom locally, whilst leaving the position and memory 
qubits unchanged. The choice of coin matrix is a function of the local 
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Figure 2 | Final variance (t = 7) of the self-avoiding quantum walk with symmetrised input state (Eq. 8), against memory recording strength 
(6 M ) and coin back-action strength (6 B ), and, as insets, some of the corresponding quantum walks (showing the full time evolution of the variance and 
probability distribution), (e) In the limit of 6 M = n/2 and 6 B = n/4 we observe ideal classical random walk statistics, since the memory is maximally 
entangled with the walker, yielding a balanced n/4 walk with decoherence upon tracing out the memory subsystem, (c) In the opposing limit of 6 M = 0 and 
irrespective of 6 B we observe ideal quantum walk statistics, since now the memory does not couple with the walker, yielding a uniform n/4 walk without 
decoherence. The memory never changes, remaining at its initial value, (a) Variance is maximised when 6 M = n/2 and 6 B = 0, for which the coin will 
always be the identity coin, yielding two outward straight line trajectories, (d) When 6 M = ^/2 and 6 B = n/2 the dispersion is restricted since the walker is 
switching between a balanced n/4 coin and a bit-flip coin at alternate steps. 



memory qubit - when q x = 0 we let =R(0c) = R(n/4) (a usual 
balanced walk with equal amplitudes of stepping in either direction), 
and when q x =lwe let U M =R(0 B ). Here, R(6) is a general PauliX 
rotation, 



R(B) = 



cos 0 



? sin 6 



sin 0 cos 6 



(4) 



Thus, 0 = 0 implements an identity gate, 6 = n/2 implements a full 
bit-flip, and 9 = n/4 implements a balanced step. Our memory will 
have the effect that if the xth qubit is set (indicating that we prev- 
iously visited location x) then the coin operator implements a par- 
tially reflecting mirror on the walker's dynamics, with strength 
determined by 0 B . This mirror is located at spatial location x and 
the reflectivity (r) is determined by 0 B (when 0 B = 0, r = 0, and when 
9 B = n/2, r = 1). Thus, the parameter 0 B controls how strong the 
effect of memory on the coin is. The probability of being reflected at a 
site where the memory is set is given by r = |sin 0 B \ 2 , otherwise it is 
given by r = |sin 9 C \ 2 . 

Propagating the walker — The step operator is defined in the usual 
way - it updates the position as a function of the coin parameter, 



S\x,c,q u ...,q N )\-->\x + c,c,qi 



,qN/ 



(5) 



The step operator affects only the position, leaving the coin and 
memory systems unchanged. 

Combined evolution — The total evolution now progresses as, 



| 1 A(f))=[S-C(0 c ,0 B )-M(0 M )]V(O)), 



(6) 



where a single step begins with memory recording, followed by the 
coin operator, and finally the step operator. 

Our formalism differs from that presented in Ref. 29, where no 
memory was required to ensure self- avoidance. Rather, the coin 
matrix was chosen so as to guarantee there is always zero ampli- 
tude in the reverse direction. This is the so-called 'non- reversing 
quantum walk'. However, such a formalism is constrained by the 
fact that one is unable to tune between different memory strengths, 
making our formalism much more versatile, as we will see in the 
next section. 

Whilst the presented formalism is in the context of a one- 
dimensional walk, it could be trivially generalised to higher dimen- 
sions, by expanding the position space. However, we will not perform 
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Figure 3 | Evolution of the self-avoiding quantum walk for 0 M = nil, 0 B — 0, with symmetrised (left) and unsymmetrised (right) input states (Eqs. 8 
and 10 respectively). This illustrates the high sensitivity of the dispersion to the initial state, as the former exhibits the maximum possible variance, 
whereas the latter exhibits zero variance. 



simulations in this context as the complexity of simulations grows 
exponentially with the number of discrete position states. 

Results 

To characterise the self- avoiding walk we will again consider the 
variance of the position probability distribution against time. We 
compute the variance, where now p x (t) is obtained by summing over 
the coin and memory degrees of freedom, 

p*w=Z) Yl Kc, qu ..., qN (t)\\ (?) 

c qi,...,q N 

effectively tracing out the coin and memory subsystems. 

In the purely classical case a 2 (t) oc t, and in the quantum case <7 2 (t) 
oc f. In our quantum walk with memory we anticipate more elab- 
orate dynamics, which we will now explore and characterise. 

In Fig. 2(centre) we plot the variance of the self-avoiding quantum 
walk against varying memory recording (6 M ) and coin back-action 
(9 B ) when evolved 7 steps. In all cases the regular coin is fixed as a 
balanced coin, 0 C = n!4. Here we have employed the symmetrised 
input state, 

|^(0)) = |l x >® -5=(|l c ) + I - l c »® \0 qi , . . . ,0^), (8) 

which guarantees symmetric position probability distributions in the 
case of the standard unbiased quantum walk. 

Each point in the plot in Fig. 2(centre) corresponds to the final 
variance of an entire quantum walk with given choices of coin back- 
action and memory recording. To understand these dynamics, the 
insets illustrate the complete quantum walk evolution dynamics cor- 
responding to some individual points in the plot. 

We observe that the variance is maximal when the memory has 
maximum recording strength (0 M = nil) and the coin has identity 
back-action (0 B = 0), as shown in Fig. 2(a). Because the memory 
recording strength is maximal, the memory qubit corresponding to 
the initial location is set. Then, because the qubit is set, the identity 
coin will be applied, which propagates the walker outwards. At the 
next step the same process applies, marking the current site and 
propagating the walker outwards. Thus, the two terms in the sym- 
metrised input state shoot off in straight line paths in opposite direc- 
tions, yielding maximum quadratic variance against time. Fitting to a 
second order polynomial, 

a 2 (t)=k 0 + k l t + k 2 t\ (9) 

the second order coefficient of the aforementioned maximum vari- 
ance walk is k 2 = 1, compared to k 2 = 0.27 for the ideal quantum 
walk, reflecting the significantly faster rate of spread of the 0 M = nil, 
0 B = 0 self- avoiding quantum walk compared to the normal ideal 
quantum walk. 



Intuitively, the maximum variance in this case appears to be a 
result of the choice of symmetrised input state. We expect that if 
we instead adopt an unsymmetrised input state, 

\il,(Q)) = \l x ,l c ,0 qi ,...,0 qN ), (10) 

there will be only a single straight line trajectory in one direction, 
yielding zero variance. This is confirmed in Fig. 3. Thus, we see that 
the variance of the self- avoiding quantum walk is highly dependent 
on the initial choice of input state, and can in fact swap us between 
maximal and zero variance. Fig. 4 illustrates the variance of an 
unsymmeterised input state across all parameter regimes. 

In the intermediate case of 0 < 0 B < nil, the coin creates a 
superposition, yielding lower variance than Fig. 2(a). 

In the case of 0 M = n!4 and 0 B = n!4 we observe statistics that are 
qualitatively very classical (Fig. 2(b)). This is because now the coin is 
a n!4 rotation irrespective of the memory (recall that 0 C = n!4 also), 
yielding a balanced walk, regardless of history. But simultaneously 
the memory register is highly entangled with the position degree of 
freedom, which, upon being traced out, yields decoherence, which 
can cause a transition from a quantum to a classical walk. This 
transition was nicely experimentally demonstrated by Broome et 
aV in the optical context. 



t=l,6 c =nl4 




Figure 4 | Variance of the self-avoiding quantum walk with 
unsymmetrised input state (Eq. 10). Contrast the different dynamics with 
Fig. 2 where a symmetrised input state was employed. In particular, note 
the drastically different variance at 6 M = nil, 6 B = 0 compared to the 
symmetrised case - the maxima in the symmetrised plot corresponds to the 
minima of the unsymmetrised plot. Note the different viewing angle than 
Fig. 2. 
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If we increase the memory recording strength to 0 M = n/2, whilst 
retaining 0 B = n/4, the walker is maximally entangling with the 
memory, yielding complete decoherence upon tracing out the mem- 
ory. Thus, in this instance we observe the binomial distribution and 
linear variance associated with the ideal classical walk (Fig. 2(e)). 

When 0 M = 0 there is no memory recording and so no entangle- 
ment between the memory and the walker. The memory stays per- 
manently fixed at its initial value and thus the coin is always the 9 C = 
n/4 coin, yielding a balanced quantum walk, irrespective of 0 B 
(Fig. 2(c)). Thus, when 0 B = n/4, the memory recording strength 
can tune us between a perfect classical walk (6 M = n/2) and a perfect 
quantum walk (0 M = 0). 

In the limiting case of maximum coin back-action (6 B = n/2) and 
maximum memory recording strength (0 M = n/2), the walker resists 
outward diffusion, yielding a distribution more localised around the 
origin (Fig. 2(d)). This is because as the walker visits sites the respect- 
ive memory qubit is completely flipped, toggling between marking 
the site as visited or not visited. When visited, the walker will incur a 
full bit-flip, sending the walker back whence it came. Otherwise a 
balanced coin is implemented and the walker spreads equally in both 
directions. Thus, on alternate steps the walker either completely flips 
direction or spreads evenly in both directions, resulting in reduced 
outward diffusion compared to a normal balanced walk. 

An interesting feature of Fig. 2 (centre) is to note that a 2 (0 B , 0 C , 0 M ) 
is not monotonically related to 0 B and 9 M . For low values of 0 B we see 
that a 2 increases monotonically with 0 M , whereas for large 0 B it 
decreases with 6 M . Similarly, for large 0 M , a 2 decreases monoton- 
ically with 0 B , whereas for small 0 M it increases with 0 M . However, 
the non-monotonic nature of a 2 may be due to the small values of t 
that we are able to simulate. In the case of unsymmetrised input states 
(Fig. 4) the dynamics are even more complex. 

Evidently, different parameter regimes for 0 B , 0 C and 9 M yield 
different couplings between the different subsystems when acted 
upon by the coin, step and memory operators. The couplings are 
illustrated in Fig. 5. The coin operator will couple all three subsys- 
tems, provided that 0 C ^ 0 B . The step operator will always couple the 
coin and position subsystems, irrespective of 0 B , 0 C and 0 M . And the 
memory operator will couple the position and memory subsystems 
provided that 0 M ^ 0. 

Comparison with classical self-avoiding random walks. In the 

classical case, self- avoiding random walks can be constructed in a 
variety of ways. For example, a classical random walk could be 
completely self- avoiding, in the sense that it may never return to a 
visited site. More generally, a metric could be defined, which specifies 
the probability with which a walker may return to a previously visited 
site. In Refs. 17,18 such a metric is defined as, 

e -gni 

where q { is the probability of stepping to site i, n { is the number of 
times site i has been previously visited, g > 0 is the strength of the 
self- avoidance (g = 0 is an ordinary classical random walk,g = o° is a 
completely self- avoiding walk, and 0 < g < o° is a weakly self- 
avoiding walk), and j sums over the nearest neighbours of the 
current vertex. Thus, the probability of returning to a given site is 
weighted exponentially with the number of times the respective site 
has been visited. Using this formalism one can tune between a regular 
non-self-avoiding classical random walk, a weakly self-avoiding 
random walk, or a completely self-avoiding random walk. A 
generic characteristic of such a self-avoiding classical random walk 
is that as g increases the walker is more likely to be found further away 
from the origin. 

The self- avoidance metric given in Eq. 1 1 biases how likely the 
walker is to revisit a site given its prior visitation history and the 
strength of the self-avoidance, g. However, this bias requires the total 



Step operator 




Memory operator 

Figure 5 | Illustration of the couplings between the position, coin and 
memory subsystems, depending on a B , 0 C and 0 M . The coin operator 
couples all three subsystems, provided that 6 C ^ 0 B . The memory operator 
couples the position and memory subsystems, provided that 6 M ^ 0- The 
step operator couples the position and coin subsystems, irrespective of the 
parameters. 

number of previous visits to each site. To compare with our quantum 
self- avoiding walk we restrict to the case where 0 M = n/2. Then the 
memory is marked and then unmarked at subsequent visits, i.e. the 
bias effectively depends on the number of visits modulo 2. Therefore, 
we modify Eq. 1 1 to count modulo 2, such that, 

e -g(nj mod2) 

C\[ = ? r- (12) 

E,-e~^ mod2 ) 

Alternately, one could achieve a situation more analogous to Eq. 1 1 
by replacing the qubits in the self- avoiding quantum walk with qu-d- 
its and letting the memory recording implement q t — > (q { + 1) mod d, 
thus giving us the ability to record up to d instances of previous visits. 
However we will not explore this option given the computational 
intractability. 

As discussed earlier, in the quantum case when 0 M = n/2 we have 
maximum coupling between the walker and the memory, yielding a 
decoherence effect. Thus, we expect the comparison between the 
classical self- avoiding random walk and the quantum self- avoiding 
walk to be closest in the regime where 0 M = n/2. 

Recall that in the classical limit o 2 (t) oc t, whereas in the quantum 
limit 0 2 (t) oc t 2 . Thus, to characterize the classical self- avoiding 
random walk we fit the variance to a 2 (t) = t fi . In Fig. 6 (top) we plot 
the exponent p against the self- avoidance strength g for a classical 
self- avoiding random walk. Contrast this with Fig. 6 (bottom), where 
we plot the analogous situation for the quantum self- avoiding walk 
with 9 M = n/2 while varying 0 B , corresponding to the classical 
regime of maximum coupling with the memory. Now 0 B takes the 
role of the self- avoidance strength. In the limit of 0 B = n/4 we observe 
an ideal classical random walk with = 1, which corresponds to 
Fig. 2(e). In the opposing limit of 0 B — 0 we observe maximum 
dispersion and ft = 2, corresponding to Fig. 2(a). Thus, within this 
region we observe a correspondence between g (in the classical case) 
and 0 B (in the quantum case). Specifically, g = 0 corresponds to 6 B = 
n/4, and g^>0 corresponds to 0 B = 0. In these two limits the self- 
avoiding classical random walk exhibits dynamics close to the 
respective self-avoiding quantum walk (Figs. 2(e) and (a) respect- 
ively). In both cases /? exhibits a monotonic relationship to the self- 
avoidance strength, ranging between /? = 1 (linear spread) and = 2 
(ballistic, quadratic spread). 

Conclusion 

We have considered a one- dimensional quantum walk with tunable 
degrees of self- avoidance. Here the walker is complemented by a 
memory register that records the previous position history and the 
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Figure 6 | Exponent of a ^(f) = f fit of the variance against time for a 
self-avoiding walk. /? = 1 corresponds to linear variance against time, 
whereas = 2 corresponds to quadratic variance against time (ballistic 
spread), (top) A classical self-avoiding random walk with self- avoidance 
function given by Eq. 12, where g is the self- avoidance strength. The fitting 
is calculated over 200 time steps and 1000 repetitions. When g = 0 we 
reduce to a normal balanced classical random walk, which exhibits linear 
spread and thus /? = 1 (analogous to Fig. 2(e)). As the self- avoidance 
strength increases, we approach a perfectly self-avoiding random walk, 
exhibiting maximum quadratic dispersion, and (3 = 2 (analogous to 
Fig. 2(a)). (bottom) The same fit for a quantum self-avoiding walk where 
6 C = ft/4 and 6 M = n/2. This is the classical regime of maximum coupling 
between the walker and the memory. Now 6 B takes the role of the self- 
avoidance strength. 

coin operator is a function of the memory, allowing the walker to 
preferentially avoid previously visited sites. 

Similar behaviour could be modelled using the formalism of Ref. 
26, where the walker has memory of past coin values rather than past 
position values. In that formalism the walker's movement may be an 
arbitrary unitary function of the walker's history, allowing arbitrary 
dynamics to be implemented. 

Whilst we have focused on one-dimensional walks, the formalism 
may be easily generalized to higher dimensions. But in that instance 
numerical simulation faces the obstacle that computational complex- 
ity grows exponentially with the number of positions. 

We have defined our self-avoiding walk such that the walker's 
direction is made a function of whether the present site has been 
previously visited. A possibility we have not explored is where the 
coin is a function of whether the neighbour in the direction of move- 
ment has been visited, but such rules often are difficult to implement 
in a unitary fashion. A variation where the walker's movements 
depends on both neighbours, though not consistent with our form- 
alism, can be made unitary and constitutes a type of quantum cellular 
automata. 

With appropriate choices of conditional coins and recording 
strengths, the walk can be made to exhibit zero dispersion or max- 
imal dispersion, which is highly dependent on the symmetrisation of 
the input state. The walker can be made to reproduce ideal quantum 



or classical statistics, or to exhibit richer, more diverse diffusive char- 
acteristics, which can be manipulated with the strength of the mem- 
ory recording (6 M ) and back-action (0 B ). In the regime of maximum 
memory recording strength, we observe a close correspondence 
between the quantum and classical self-avoiding walks, where the 
coin back-action takes the role of the self- avoidance strength. 

The rich diffusive characteristics of the quantum self- avoiding 
walk may be applicable to modelling a broader range of diffusive 
phenomenon than ordinary quantum walks, and may lend them- 
selves to novel quantum information processing applications such as 
applications in quantum simulation. 
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